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ABSTRACT 

Gravitational microlensing at cosmological distances is potentially a powerful tool 
for probing the mass functions of stars and compact objects in other galaxies. In the 
case of multiply-imaged quasars, microlensing data has been used to determine the 
average microlens mass. However the measurements have relied on an assumed trans- 
verse velocity for the lensing galaxy. Since the measured mass scales with the square 
of the transverse velocity, published mass limits are quite uncertain. In the case of 
Q2237-I-0305 we have properly constrained this uncertainty. The distribution of light 
curve derivatives allows quantitative treatment of the relative rates of microlensing 
due to proper motions of microlenses, the orbital stream motion of microlenses and 
the bulk galactic transverse velocity. By demanding that the microlensing rate due to 
the motions of microlenses is the minimum that should be observed we determine lower 
limits for the average mass of stars and compact objects in the bulge of Q2237-I-0305. 
If microlenses are assumed to move in an orbital stream the lower limit ranges be- 
tween 0.005 and 0.023AfQ where the the systematic dependence is due to the fraction 
of smooth matter and the size of photometric error assumed for published monitoring 
data. However, if the microlenses are assumed to move according to an isotropic ve- 
locity dispersion then a larger lower limit of 0.019-0. IIM© is obtained. A significant 
contribution of Jupiter mass compact objects to the mass distribution of the galactic 
bulge of Q2237-I-0305 is therefore unambiguously ruled out. 

Key words: gravitational lensing - microlensing - stellar masses. 



1 INTRODUCTION 

Q2237+0305 comprises a source quasar with a redshift of 

2 = 1.695 that is gravitationally lensed by a foreground 
galaxy at z — 0.0394 producing 4 images with separations 
of ~ 1". Each of the 4 images are observed through the 
bulge of a galaxy which has an optical depth in stars that is 
of order unity (eg. Kent & Falco 1988; Schneider et al. 1988; 
Schmidt, Webster & Lewis 1998). In addition, the proxim- 
ity of the lensing galaxy means that the effective transverse 
velocity may be high. The combination of these facts make 
Q2237-I-0305 the ideal object from which to study microlens- 
ing. Indeed, Q2237-I-0305 is the only object in which cosmo- 
logical microlensing has been confirmed (Irwin et al. 1989; 
Corrigan et al. 1991). 

Numerical microlensing simulations (eg. Wambsganss, 
Paczynski & Schneider 1990; Witt, Kayser & Refsdal 1993) 
have shown that the statistics of microlensed high magnifi- 
cation events (HMEs) obtained from long term monitoring 
may provide information on properties of the lens such as 
the stellar mass function and the percentage of mass in stars 



for the bulge. However, the monitoring period required is 
greater than 100 years (Wambsganss, Paczynski & Schnei- 
der 1990). There are several other unknown quantities in the 
problem. These include the magnitude and direction of any 
transverse motion, as well as the source size. In particular, 
the magnitude and direction of the transverse motion are 
degenerate with the density of caustics (a function of the 
mean compact object mass) for a given set of microlensing 
statistics. 

Microlensed fiuctuation in the quasars continuum re- 
sults from motion due to both a galactic transverse velocity 
and to the random proper motion or stream motion of stars 
and compact objects. Wyithe, Webster & Turner (1999a) 
(hereafter WWTa) define the equivalent transverse veloc- 
ity as the transverse velocity in a model containing stars 
with static positions, that produces a rate of microlensing 
most closely resembling that of a model containing stellar 
proper motions. Foltz et al. (1992) measured the central ve- 
locity dispersion of Q2237-I-0305 to be ~ 215km sec~^ , and 
theoretical models (Schmidt, Webster & Lewis 1998) pre- 
dict a value of ~ 165A;m sec~^ . If the dispersion is isotropic 
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then (as discussed in Wyithe, Webster & turner (1999b) 
(hereafter WWTb)) the equivalent transverse velocity cal- 
culated from the overall microlensing rate (all 4 images) of 
Q2237+0305 is larger than this value. The microlensing rate 
that results from the line-of-sight velocity dispersion (of an 
isotropic distribution) is therefore comparable to that of the 
likely transverse velocity. This means that the often made 
assumption that random proper motions provide a negligible 
contribution to microlensing in Q2237-I-0305 is incorrect. 

While the microlensing rate resulting from random stel- 
lar motions is dependent on the mean microlens mass, the 
equivalent transverse velocity (WWTa; WWTb) that de- 
scribes the microlensing rate is not. We use this fact to 
break the degeneracy between microlensing rate and mean 
microlens mass, and hence obtain useful limits on the mass 
function along the line-of-sight through the bulge of the lens- 
ing galaxy. 

This paper is presented in 5 parts. Sec. H contains a 
general discussion of microlensing and the mass function. 
Sec. ^ discusses the numerical methods used to model mi- 
crolensing in Q2237-I-0305 and Sec. ^ describes a method to 
place limits on the mass function through consideration of 
microlensing due to both a transverse velocity and a stellar 
velocity dispersion. 



2 MICROLENSING AND THE STELLAR 
MASS FUNCTION 

Stellar mass functions have traditionally been measured 
through the combination of an observed luminosity function 
and an empirical mass-luminosity relationship. However, be- 
cause of the difficulties inherent in the observations of faint 
stars, these determinations become uncertain as the hydro- 
gen burning limit (~ O.OSM©) is approached and crossed. In 
contrast to this approach, microlensing uses the mass of a 
lens to magnify the observed flux of a background source 
through gravitational deflection of the light bundle. The 
statistics obtained are therefore free of any bias introduced 
by the hydrogen burning limit, and so microlensing is a pow- 
erful tool for determining the contribution to the mass func- 
tion of low mass stars and dark compact objects. 

Gravitational microlensing is observed in two very dif- 
ferent regimes. Paczynski (1986) suggested surveying mil- 
lions of stars in the Magellanic Clouds for microlensing in- 
duced flux variation as a way of detecting solar mass com- 
pact objects in the halo of our own galaxy. Several groups 
have undertaken such searches (eg. Alcock et al. 1997a; Re- 
nault et al. 1998). In a similar vein, Paczynski (1991) and 
Griest et al. (1991) suggested microlensing experiments to- 
wards the Galactic bulge as a method to probe the masses of 
stars along the line-of-sight in the Galactic disc. Microlens- 
ing searches towards the Galactic bulge have since found a 
higher rate of microlensing events than are found towards 
the LMC (eg. Alcock et al. 1997b; Udalski et al. 1994). In 
a very different microlensing regime, the observed flux of a 
back-ground quasar can be altered by gravitational lensing 
by stars or halo objects in a foreground galaxy. Q2237H-0305 
is an example of a quasar that lies very close to a galactic 
line-of-sight. In such cases multiple imaging occurs, allowing 
microlensed variation to be easily separated from intrinsic 
fluctuations which are observed in all images. 



The analyses of microlensing in the Galactic and cosmo- 
logical regimes are very different. While there are only a few 
sources (quasar images) in the cosmological case rather than 
the millions available in Galactic microlensing searches, the 
optical depth (or equivalently the probability of lensing) is 
~ 10® greater. The transverse velocity is a critical param- 
eter for the determination of a mass or masses responsible 
for a microlensing event, regardless of the microlensing sce- 
nario. Unfortunately it is unknown in both cases. In Galactic 
microlensing calculations the velocities are inferred from an 
assumed distribution. However in the cosmological case the 
transverse velocities of the lensing objects are not indepen- 
dent. Rather, they are equal except for the contribution of 
the individual stellar proper motions. Moreover, cosmologi- 
cal microlensing results from the gravitational contribution 
of a large ensemble of many hundreds of masses rather than 
on a single microlens. Microlensed, multiply-imaged quasars 
are an excellent probe of the mass function of compact ob- 
jects along the quasar image line-of-sight because they in- 
teract with a large number of microlenses. 

Several authors have placed limits on the masses of mi- 
crolenses responsible for cosmological microlensing. Schmidt 
& Wambsganss (1998) use the lack of observed variation in 
Q0957+561 to place a lower limit on the mass of microlenses 
in the halo of the lensing galaxy. They rule out a mean mass 
of halo objects (m) <C 1Q~^Mq. However their determina- 
tion is dependent on the source size, and the fraction of 
smooth matter employed in their calculations. A transverse 
velocity of Vt = 600 fern sec ~^ is assumed. The uncertainty 
in this assumed value is the most serious problem for this 
kind of analysis because the values of mass obtained are 
cx «t . Lewis & Irwin (1996) compared the monitoring data 
of Q22374-0305 with simulations using a structure function 
to analyse variability. They too assume a transverse veloc- 
ity of Vt = 600 km sec ~^ and conclude that the mean mass 
of objects in Q2237+0305 is O.IM© ^ (m) <, IOMq. This 
result is also proportional to the square of the unknown 
transverse velocity. 

While they do not measure precisely the same quantity, 
the above results are supportive of those of the MACHO 
experiment (eg Alcock et al. 1997a,b). From microlensing 
observations towards the LMC by the Galactic halo, Alcock 
et al. (1997a) conclude that the average MACHO masses are 
0.08 — 0.93Mq. The quoted range includes both statistical 
uncertainties and systematic effects introduced by the as- 
sumption of halo model. In addition, they flnd evidence for 
an excess of events over those predicted from stellar lensing 
alone. Alcock et al. (1997b) flnd that the range of time-scales 
towards the Galactic bulge are consistent with microlens 
masses of O.IMq ^ (m) ^ I.OMq. Tab. |l| summarises the 
results described. 



3 THE MICROLENSING MODELS 

3.1 Microlensing models for Q2237+0305 

Throughout the paper, standard notation for gravitational 
lensing is used. The Einstein radius of a IMq star in the 
source plane is denoted by rjo. The normalised shear is de- 
noted by 7, and the convergence or optical depth by k. The 
model for gravitational microlensing consists of a very large 
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Table 1. A collection of results from microlensing data for the average mass of compact objects in galactic hales and/or bulges of the 
Milky Way and Quasar lensing galaxies. 



Milky Way halo 
(Alcock et al. 1997a) 


Results for mean masses of dark halo objects toward LMC. 
The ranges include statistical (lo") as well as systematic 
uncertainties due to choice of halo model. 


(m) = O.lSl 
-0.551 




A/Till , "WT^ l-,.,!™^ 

Miiky Way bulge 
(Alcock et al. 1997b) 


Microlensing event time-scales are consistent with a mean 
mass of compact objects and stars in the Milky Way bulge. 


O.IMq <, (m> <, 1.OM0 


0957+561 halo 
(bchmidt &; Wambsganss 1998) 


The limit has 99% and 95% significance levels for assumed 
source sizes oi 4 X 10 cm and 4 X 10 cm, and scales as 
the square of the assumed transverse velocity vt- 


(m> ^ 0.001 


\mo 1 


2237+0305 bulge and halo 
(Lewis & Irwin 1996) 


The microlensing rate of current light curves is consistent 
with a mean mass of stars and compact objects. The value 
scales as the square of the assumed transverse velocity vt- 


O.lM0^(m>(^)'<:iOM0 


2237+0305 bulge and halo 
This paper. 


The proper motions of microlenses provide a minimum micro- 
lensing rate. The limits given for the mean mass have a 99% 
significance level, and depend on the bulge dynamics assumed. 


(m) = O.OlOl 
-0.291 





Table 2. Values of the total optical depth and the magnitude of 
the shear at the position of each of the 4 images of Q2237+0305. 
The quoted values are those of Schmidt, Webster & Lewis (1998). 
K» and Kc are the optical depths in stars and in smoothly dis- 
tributed matter respectively. 



Image 




|7| 


A 


0.36 


0.40 


B 


0.36 


0.40 


C 


0.69 


0.71 


D 


0.59 


0.61 



sheet of point masses that simulates the section of galaxy 
along the image line-of-sight, together with a shear term 
that includes the perturbing effect of the mass distribution 
of the lensing galaxy as a whole. The normalised lens equa- 
tion for a field of point masses with an applied shear in terms 
of these quantities is 

Here x and y are the normalised image and source posi- 
tions respectively, and the and m-' are the normalised 
positions and masses of the microlenses. Kc is the optical 
depth in smoothly distributed matter. Eqn. ^ is solved for 
the macroimage magnification at many points along a pre- 
defined source trajectory through the inversion technique of 
Lewis et al. (1993) and Witt (1993). The region of the lens 
plane in which image solutions need to be found to ensure 
that 99% of the total macro-image flux is recovered from 
a source point was described by Katz, Balbus & Paczynski 
(1986). The union of areas of in the lens plane corresponding 
to the flux collection area of each point on the source line 
is known as the shooting region, the method for determin- 
ing the dimensions of which is described in Lewis & Irwin 
(1995), and Wyithe & Webster (1999). The radius of the 



disc of point masses is chosen to be 1.2 times that required 
to cover this shooting region. 

We assume the macro-parameters for Q2237+0305 cal- 
culated by Schmidt, Webster & Lewis (1998) (the values are 
shown in Tab. Two models are considered for the mass 
distribution, one with no continuously distributed matter 
and one where smooth matter contributes 50% of the surface 
mass density. Both the microlensing rate due to a transverse 
velocity (Witt, Kayser & Refsdal 1993; Lewis & Irwin 1996; 
WWTb), as well as the corresponding rate due to proper 
motions (WWTa) are not functions of the detail of the mi- 
crolens mass distribution, but depend only on the mean mi- 
crolens mass. The mean of the microlens mass function can 
therefore be determined independently from its form. On the 
other hand, information on the form of the mass function 
is unobtainable through consideration of the microlensing 
rate. We limit our attention to models in which all the point 
masses have identical mass. 

Each of our models was computed for a source track 
of length lOrjo (corresponding to ~ 70 years for an effec- 
tive transverse velocity of 600 km sec~^). We set 500 to be 
the minimum number of stars in a model. Two orienta- 
tions were chosen for the transverse velocity with respect 
to the galaxy, with the source trajectory being parallel to 
the A— B or C— D axes. At each orientation, 100 simulations 
were made for each of the 4 images of Q2237+0305 in com- 
bination with each of the model mass distributions. The two 
orientations bracket the range of possibilities, and because 
the images are positioned approximately orthogonally with 
respect to the galactic centre correspond to shear values of 
7A,7s < 0,7c, 7d > and 7a,7s < 0,7c,7d > respec- 
tively. Tab. ^ shows the number of stars used for each of 
these models along with the mean magnification (/j.) , and the 
theoretical magnification {fj,th) for comparison. The average 
magnification was found from the combination of the two 
sets of models, and the error calculated from the standard 
deviation of a subset of 6 simulations (3 per model orienta- 
tion) . Fig. ^ shows the magnification distributions for each 
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Table 3. Values for the number of stars in the microlensing mod- 
els (left column 7 > 0; right column 7 < 0), the average model 
magnifications and the theoretical values for comparison. 



0% smooth matter 



Image 


No. Stars 






A/B 
C 
D 


763 500 
1022 500 
2887 1301 


3.97±.28 
2.39±.ll 
4.79±.14 


4.01 
2.45 
4.90 


50% smooth matter 


Image 


No. Stars 




ifJ-th) 


A/B 
C 
D 


500 500 
500 500 
1080 500 


4.01±.19 
2.44±.09 
4.84±.20 


4.01 
2.45 
4.90 
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Figure 1. The magnification distributions for images A/B (left), 
C (centre) and D (right). The top and bottom rows show the 
distributions for models that have no smooth matter component, 
and a 50% smooth matter component respectively. The light and 
dark lines represent distributions that are sampled parallel to and 
at right angles to the shear. 

image in each model. Simulations having 7 < and 7 > 
are represented by hght and dark lines respectively. The dis- 
tributions are quantitatively comparable to those presented 
for similar values of ft and 7 in Lewis & Irwin (1995), and 
demonstrate that the magnification distribution is indepen- 
dent of the source direction as required. 

We assume that microlensing is produced through the 
combination of a galactic transverse velocity with each of 
two classes of proper motion for individual stars with re- 
spect to the galaxy: an isotropic velocity dispersion and a 
circular stream motion. Moreover, we assume that the mag- 
nitude of the dispersion or stream motions are the same for 
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Effective Transverse Velocity (km/sec) 

Figure 2. Plots of the cumulative probability (Pv) for the ef- 
fective transverse velocity. The solid and dot-dashed lines corre- 
spond to the models with trajectory directions that are aligned 
with the C-D and A-B image axes respectively. The dark and 
light lines represent results from models containing 0% and 50% 
continuously distributed matter. The upper and lower plots cor- 
respond to the cases where the photometric error was assumed to 
be (T5£;=0.01 mag in images A/B and o'5£=0.02 in images C/D, 
and o'i^E=0.02 mag in images A/B and o'i^E=0.04 in images C/D. 
The model assumed \Mq microlenses. 

each of the four images. We take the theoretical value of 
(T, ~ 165 fem see" ^ for the line-of-sight velocity dispersion 
of the stars in the galactic bulge. This is lower than the 
value observed by Foltz et al. (1992), so lower limits placed 
on the mass are more conservative. 

3.2 The effective transverse velocity 

We define the effective transverse velocity as the transverse 
velocity that produces a microlensing rate from a static field 
model equal to that of the observed light-curve. The effective 
transverse velocity therefore describes the microlensing rate 
due to the combination of the effects of a galactic transverse 
velocity and random microlens proper motion. 

Our calculation of effective transverse velocity utilises 
the cumulative histogram of derivatives calculated from the 
6 difference light-curves (A-B, A-C, A-D, B-C, B-D, C-D). 
Through computation of the difference light-curves, intrin- 
sic flux variation, as well as the systematic component of 
the observational uncertainty are removed from the data. 
WWTb describes a procedure for determining the proba- 
bility that the effective transverse velocity is less than an 
assumed value. That paper provides correct upper limits, 
but does not obtain the cumulative probability function. A 
modifled procedure to determine the cumulative probabil- 
ity for effective transverse velocity from ensembles of mock 
data-sets is described below. 

Many mock observations of the 6 difference light-curves 
were produced from model light-curves with a sampling rate 
and period that are identical to the published monitoring 
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data (Irwin et al. 1989, Corrigan et al. 1991, 0stensen et al. 
1996). Observational errors were also simulated by apply- 
ing a random fluctuation to each point on the light curve, 
distributed according to a Gaussian with a half-width a de- 
scribing the published photometric error. The simulations 
used two different estimates for the uncertainty in the pho- 
tometric magnitudes. In the first case a small error was as- 
sumed (SE). For images A and B, asE=0.01 mag, and for 
images C and D asE=0.02 mag. In the second case, a larger 
error was assumed (LE). For images A and B, (7l_e=0.02 
mag and for images C and D aLE—O.Oi mag. The observa- 
tional error in Irwin et al. (1989) was 0.02 mag. We note 
that the adoption of a larger error leads to a smaller bound 
on the measured upper limit of effective transverse velocity, 
and subsequently a larger lower bound on the mean mass. 

The simulations assume a point source since the value of 
effective transverse velocity measured with a model is not a 
sensitive function of the source size (in the range of interest) 
assumed (WWTb). The insensitivity arises from the fact 
that the measurement is derived predominantly from points 
that are not part of a HME. 

We define a set of effective transverse velocities (Ve//) 
and determine the probability, based on light-curve deriva- 
tives that each of these correctly describes the true value 
for Q2237-I-0305 (Vgai)- For each model, 5000 mock obser- 
vations were produced at each effective transverse velocity 
from the four sets (one per image) of 100 IO770 simulated 
light-curves. Due to the finite sampling period, histograms 
of microlensed light-curve derivatives produced from the in- 
dividual mock observations describe typical but not average 
behaviour. Therefore an average histogram was also com- 
puted from each set of mock observations. For each mock 
observation at each pre-defined effective transverse velocity 
(Ve//), a mock measurement of effective transverse velocity 
Vef f was made by minimising the KS difference between the 
average histogram for Veff and the histogram of the mock 
observation. Thus we calculate the function representing the 
likelihood pihoodivef f\Vef f) for observing Vefj given an as- 
sumption for the true value (14//). 

Similarly, from the observed histogram of microlensed 
light-curve derivatives we find the effective transverse ve- 
locity {vobs) that best describes the observed microlensing 
rate. Using Bayes' theorem we calculate the posterior prob- 
ability that the effective galactic transverse velocity Vgai is 
less than an assumed value Veff- 

Pv{Vgal < Veff\Vobs) = 

Plhood{Vobs\Vif f) Pprtor{Vif f) dV^f f , (2) 

where Pprior{Vef f) is the prior probability for Ve//. We have 
assumed two different and physically un-motivated priors. 
These are respectively flat: 



Pprior{Veff) OC dVe// 
PprioriVeff) = 

and logarithmic: 



Pprior (Veff) 
Pprior (Veff) 



dVe 



V 







< Veff < 2000V ("i) 
otherwise, (3) 



< Ve// < 2000 V (m) 
otherwise. (4) 



These priors bracket the range and spread of physically plau- 



sible priors such as a Gaussian based on observed peculiar 
velocities (eg Mould et al. 1993) or those computed from 
numerical studies. Our results are insensitive to the choice 
of these priors, and we conclude that our estimate of proba- 
bility is dominated by the microlensing observations rather 
than our assumed prior for Ve//. 

Pv{vgai < Veff\vobs) IS plotted iu Fig. 1^ for the models 
discussed in this paper ((m) — IMq). In these plots the solid 
and dot-dashed lines correspond to source trajectories along 
the C-D image axis and A-B image axis, and the dark and 
light lines correspond to models with 0% and 50% smoothly 
distributed matter. The upper and lower plots correspond 
to the cases where the photometric error was assumed to 
be asE~0.01 mag in images A/B and as_B=0.02 in images 
C/D, and a_L_E=0.02 mag in images A/B and aLE~0.04 in 
images C/D. 



3.3 The contribution to microlensing of stellar 
proper motions 

The average histogram of light-curve derivatives can be com- 
puted for the case where microlensing results from a random 
velocity dispersion rather than a transverse velocity. We de- 
fine the equivalent transverse velocity as the transverse ve- 
locity that in combination with a static microlensing model 
produces a microlensing rate closest to the rate produced by 
a random velocity dispersion of the point masses. The mi- 
crolensing rates are considered equivalent at the transverse 
velocity that produces a cumulative histogram of light-curve 
derivatives closest (has the minimum possible KS difference) 
to the corresponding proper motion histogram. 

The formalism required for the computation of the cu- 
mulative distribution of derivatives resulting from proper 
motion of stars was developed in WWTa. The upper panel 
of Tab. ^ shows the values of equivalent transverse veloc- 
ity of a Gaussian 1-d velocity dispersion with a half-width 
of cr« = 165 km sec~^. These values were computed from 
the difference light-curves for models of Q2237-I-0305 (the 
quoted error describes the range of values obtained over 
three separate sets of simulations). The two cases shown 
are for a trajectory that is parallel to the shear in images A 
and B (7a, 7s > 0), and for one that is parallel to the shear 
in images C and D {"/c,Jd > 0). In all cases the effective 
transverse velocity is larger than the 1-d velocity dispersion. 



4 THE MASS OF COMPACT OBJECTS IN 
Q2237+0305 

4.1 Placing limits on the mass using effective 
transverse velocity 

In this section we describe a method for determining the 
lower limit for the average compact object mass. The mi- 
crolensing rate has a minimum possible value determined 
by the size of the stellar velocity dispersion, thus there is a 
minimum mass which can explain the observed rate. This is 
quantified by noting that the measured effective transverse 
velocity must be greater than the equivalent transverse ve- 
locity. 

There is a simple scaling that relates models of the same 
optical depth but consisting of model stars with a different 
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Table 4. Table showing the equivalent transverse velocities and the 99%, 95% and 90% lower limits for the mean microlens mass as well 
as the most likely value. Values for the mass limits are shown corresponding to the cases where the error was cTg^ = 0.01 mags in images 
A/B, asE = 0.02 mags in images C/D, and where the error was (Ti^e = 0.02 mags in images A/B, ctlb = 0.04 mags in images C/D (in 
parentheses). Values are shown that correspond to an isotropic velocity dispersion of a* = 165 fcmsec^^ (top table), a circular stream 
velocity of Vstream = 233 km sec~^ (bottom table). 



(J* = 165 km sec 



Trajectory 




Smooth 


Equiv. 


Bayesian 


mioyj{99%) 


m;oti,{95%) 


m;o„(90%) 


mode 


Orientation 




Matter 


Vel. (km sec^^) 


prior 


(Mq) 


(Mq) 


(Mq) 


(Mq) 


Shear: 




0% 


280 ±15 


log 


0.050 (0.111) 


0.111 (0.277) 


0.176 (0.484) 


0.139 (0.288) 


7A.7S > 0, -fc^'TD 


< 


50% 


199 ±12 


log 


0.019 (0.039) 


0.042 (0.095) 


0.066 (0.161) 


0.054 (0.087) 


Shear: 




0% 


270 ±10 


log 


0.053 (0.111) 


0.111 (0.270) 


0.173 (0.463) 


0.139 (0.281) 


7A,7s < 0, 7c,7D 


> 


50% 


181 ±4 


log 


0.018 (0.036) 


0.039 (0.084) 


0.061 (0.143) 


0.042 (0.087) 


Shear: 




0% 


280 ±15 


flat 


0.098 (0.230) 


0.270 (0.739) 


0.505 (1.496) 


0.222 (0.450) 


7A,7s > 0, 7c,7D 


< 


50% 


199 ±12 


flat 


0.042 (0.090) 


0.121 (0.301) 


0.240 (0.643) 


0.069 (0.139) 


Shear: 




0% 


270 ±10 


flat 


0.107 (0.244) 


0.293 (0.792) 


0.553 (1.638) 


0.222 (0.450) 


7A,7s < 0, 7Ci7d 


> 


50% 


181 ±4 


flat 


0.044 (0.093) 


0.135 (0.330) 


0.278 (0.729) 


0.069 (0.239) 



'(^stream — 233 fc?Tl SGC 



Trajectory 




Smooth 


Equiv. 


Bayesian 


m.io™(99%) 


m;ou,(95%) 


mio™(90%) 


mode 


Orientation 




Matter 


Vol. (kmsec^^) 


prior 


(Mq) 


(Mq) 


(Mq) 


(Me) 


Shear: 




0% 


129 ±7 


log 


0.010 (0.023) 


0.024 (0.061) 


0.039 (0.113) 


0.027 (0.043) 


7A.7fl > 0, 7Ci7d 


< 


50% 


99 ±5 


log 


0.006 (0.011) 


0.013 (0.030) 


0.022 (0.056) 


0.013 (0.021) 


Shear: 




0% 


120 ±2 


log 


(0.011 (0.023) 


0.025 (0.064) 


0.042 (0.120) 


0.021 (0.043) 


7A,7s < 0, 7C,7-D 


> 


50% 


85 ±3 


log 


0.005 (0.009) 


0.011 (0.024) 


0.018 (0.044) 


0.010 (0.021) 


Shear: 




0% 


129 ±7 


flat 


0.028 (0.071) 


0.105 (0.319) 


0.251 (0.823) 


0.034 (0.069) 


7a,7b > 0. IClD 


< 


50% 


99 ±5 


flat 


0.014 (0.031) 


0.053 (0.144) 


0.136 (0.395) 


0.017 (0.027) 


Shear: 




0% 


120 ±2 


flat 


0.027 (0.069) 


0.108 (0.33) 


0.273 (0.897) 


0.034 (0.687) 


1A,1B < 0, 7C,7-D 


> 


50% 


85 ±3 


flat 


0.014 (0.031) 


0.060 (0.159) 


0.162 (0.454) 


0.017 (0.027) 



mass. The dimensionless Einstein radius of a point mass is 
y/m. Therefore, reducing all microlens masses by a factor a 
reduces all physical distances (in both the source and lens 
planes) between dimensionless coordinates by a factor ^/a 
(eg. Witt, Kayser & Refsdal 1993). For a given transverse 
velocity this results in an increase by a factor of 1 in the 
gradient of the light-curve at all points. Similarly, the rate 
of change of magnification at a static source point due to 
stellar proper motions increases by a factor Xj^/a. Since it 
results from the proper motions the minimum rate is larger 
if a smaller mean microlens mass is assumed. However the 
minimum rate cannot statistically exceed the observed rate. 
The ratio of equivalent transverse velocity to velocity disper- 
sion is independent of the mean microlens mass, however the 
measurement of effective transverse velocity is proportional 
to ^/a. This allows limits to be placed on the value of a and 
therefore the stellar mass from the observed microlensing 
rate. 



4.2 Mass limits from random proper motions 

For a model where 7a, 7s > 0, and the simulated errors are 
assumed to be cts_b = 0.01 and asE = 0.02 mags in images 
A/B and C/D respectively, the effective transverse velocity 
is less than 1346y^ (m) kmsec~^ at the 99% level, where (m) 
is the microlens mass assumed. In the absence of a galac- 
tic transverse velocity, the minimum level of microlensing is 
described by the equivalent transverse velocity of the stel- 
lar velocity dispersion Vequiv{DISP) = 280 fcm see" ^ (see 
Tab. ^). The value of (m) at which Kipper (99%) drops be- 
low Vegmv{DISP) IS m,<,„(99%): 

The minimum mass of stars in this model is therefore ~ 
0.04Mq . This calculation assumes a value for the one dimen- 
sional velocity dispersion of 165 km sec~^ . However we note 
that all lower limits obtained in this section are proportional 
to the square of this value. The procedure also assumes that 
the measured effective transverse velocity is independent of 
the assumed source size. WWTb show that this is approxi- 
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Mean Microlens Mass (solar masses) Mean Microlens Mass (solar nnasses) 

Figure 4. Contours of percentage peak height of the function pm.v{{nfi) , Vtran)- The contours shown are the 0.1%, 1.0%, 3.6%, 14%, 26% 
and 61% levels. The solid and dotted curves represent the resulting functions when the photometric error was assumed to be a-g£;=0.01 
mag in images A/B and (Tg£=0.02 in images C/D, and (Ti£;=0.02 mag in images A/B and (Ti£;=0.04 in images C/D. 



mately true for source sizes smaller than « 4 x 10 ymcm. 
This value compares favourably with the findings of Wamb- 
sganss, Paczynski & Schneider (1990) that a source smaller 
than ~ 2 X 10^^ ^(m)/. 225 cm is required to typically re- 
produce the observed HME amplitude. 

The curves shown in Fig. Q represent the cumulative 
probability PviVgai < Veff\/rn) that the effective transverse 
velocity Vgai is less than Vgffy'Tn. It follows that the cumu- 
lative probability that the mean microlens mass is less than 
the value 



M = 



Veff 



Um{l,m) < M\V,gu^v) = 1 ~ PviVgal < Vef fV M) 



(6) 



(7) 



The calculation of mjo^ (99%) performed above assumes 
that the physical transverse velocity of the galaxy with re- 
spect to the observer source line of sight Vtran = 0. How- 
ever since the galaxy may have a component of transverse 
motion, Vequiv, which describes the minimum microlensing 
rate is a function of Vtran- The dependence of Vequiv on 
Vtran is determined according to the relationships described 



in WWTb. Since the galactic transverse velocity is not nec- 



essarily zero, Um{{m) < M\Vet^ 



must be determined at 



a series of transverse velocities and combined with prior 
probabilities for Vtran to obtain estimates for the proba- 
bility of the mean microlens mass. We have assumed the 
flat and logarithmic priors employed for the calculation of 

P.iVgal < Veff): 



Um{{m) < M) 



Um.{{m) < M\Vtran)Pprior{Vtran)dVtrar 



(8) 



The probability density pm{{m)) for the mean microlens 



mass is obtained by taking the derivative 



d (m) 



Note that 



this distribution does not represent the mass function. 

Fig. ^ shows the functions pm{{m)) for the 0% smooth 
matter and 50% smooth matter models for both directions 
considered. The solid and dotted curves represent functions 
assuming the photometric error to be asE=0.01 mag in im- 
ages A/B and crs_B=0.02 in images C/D, and aLE~0.02 
mag in images A/B and (t_l_b=0.04 in images C/D. The 
light and dark lines refer to the logarithmic and flat priors 
for transverse velocity. Results for rr!,;ou,(99%), (95%), 
miom{90%) and the mode of the distribution are shown in 
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Mean mass (solar masses) 

Figure 3. The difTcrontial distributions for tlie mean microlens 
mass. Tlie solid and dotted curves represent the resulting func- 
tions when the photometric error was assumed to be crs^=0.f)l 
mag in images A/B and CT5e=0.02 in images C/D, and cr/^£;=0.02 
mag in images A/B and (T££;=0.04 in images C/D. The light and 
dark lines correspond to the assumptions of logarithmic and flat 
priors for effective transverse velocity referred to in the text. 



the top panel of Tab. |^ In each column, the values cor- 
respond to the two assumptions for the simulated error 
referred to above. m;otu (99%) is > O.O2M0 in all models 
considered. The mode of the distribution ranges between 
O.O4M0 and O.45A/0. It is interesting to note that these 
models produce distribution modes that are similar to the 
mean of a Salpeter distribution with a lower cutoff at the 
hydrogen burning limit. We also note that the distribution 
is highly asymmetric, with the mode being of the same order 
as the lower limit. 



4.3 The effect of systematic assumptions 

In contrast to the calculation of Pv{Vgai < Veff \/ {m)) the 
assumed prior for Ve// makes a significant difference to the 
microlens mass limits obtained. In particular, low mass mi- 
crolenses are less likely if large transverse velocities are as- 
sumed. This dependence on the prior illustrates the limi- 
tation of the light curve data for breaking the degeneracy 
between galactic transverse velocity and microlens mass in 
the high mass- velocity regime. 

The assumption of a larger photometric uncertainty 
yields a lower estimate of the transverse velocity (WWTb). 
This in turn means that a larger estimate is made for the 
mean microlens mass. The effect is readily apparent in Fig. ^ 
as well as from Tab. ^. 

We find lower microlens mass limits if a non-zero frac- 
tion of smooth matter is assumed. This results from the 
combination of a larger determination of effective transverse 
velocity and a smaller relative contribution to microlensing 
of the stellar proper motions resulting in a smaller minimum 
for microlensing rate. This is expected since no microlensing 
will be observed due to a galaxy composed entirely of contin- 
uously distributed matter. The plots in Fig. ^ demonstrate 
the level of dependence on smooth matter component. 

Pm{{'m)\Vtran) Is extremely sensitive to the level of 
photometric error assumed at large (m). In addition, mi- 
crolenses of arbitrarily high mass can produce the required 
rate in combination with a correspondingly large transverse 
velocity. Therefore at large (m), pm({m)) is a sensitive func- 
tion of the prior assumed. For these reasons our method does 
not place useful upper limits on the mean microlens mass. 



4.4 Simultaneous limits on mean microlens mass 
and transverse velocity 

To explore the co-dependence of the determined mean mi- 
crolens mass and galactic transverse velocity we have com- 
puted the two dimensional distribution 

ran) — pm 

{{■m)\Vtran)Pv{Vtran\{m)). (9) 

Figure ^ shows contour plots of this distribution for the as- 
sumed source orientations, smooth matter fraction and pho- 
tometric errors. All possible combinations of the two values 
of each parameter are included, so that eight models are 
shown. The contours are at 0.1, 1.0, 3.6, 14, 26, 61 percent 
of the peak height, and so the extrema of the inner 4 contours 
represent Aa, Scr, 2a, and la limits on the single variables. 
The most important benefit of combining the distributions 
of Vtran and (m) is that upper limits can be placed on the 
mean microlens mass and galactic transverse velocity. The 
limits on (m) and Vtran are subject to both random and sys- 
tematic uncertainties. Assuming that the chosen values for 
these parameters bracket the likely models, upper and lower 
limits for (m) , and an upper limit for Vtran can be found by 
using the most extreme values of the contour limits. This 
yields the 95 per cent result that O.OIM© <, (m) <, IMq, 
and Vtran ^ 300 fcmsec"^ (95%). The upper limit is lower 
than that obtained by Lewis & Irwin (1996), although suffi- 
ciently high that it does not conflict with any popular mod- 
els of galactic halos or bulges. The elongation of the outer 
contours with a power law slope of ^ illustrates the degener- 
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effective transverse velocity increases. This effect is mag- 
nified by the inclusion of a component of smooth matter. 
In addition, the equivalent transverse velocity of the stellar 
proper motions is also increased under the same circum- 
stances. Interestingly, during our calculation of lower mass 
limit these effects tend to have a cancelling effect, render- 
ing the limits reasonably insensitive to the source trajectory 
direction assumed. 



4.6 Mass limits from stellar stream motions 



Einsteiti cadii 



Emsteln cadll 



Figure 5. The magnification maps for images A (top) and C 
(bottom), in the cases of 0% smooth matter (left) and 50% smooth 
matter (right). The point masses each had a mass of IMq 



acy between microlens mass and transverse velocity at large 
mass. 



4.5 The effect of trajectory direction 

Fig. ^ shows the magnification maps for images A and C 
in the cases of 0% and 50% continuously distributed mat- 
ter. The caustic networks in both images are significantly 
stretched so that a source could travel a significant distance 
in the direction of caustic clustering and be subject to lit- 
tle microlensed fiux variation. The measured effective trans- 
verse velocity is dependent on trajectory direction, with a 
larger measurement resulting from cases where the source 
trajectory is perpendicular to the shear in images C and D 
(7c, 7d < 0). The larger measurement quantifies the greater 
stretching of the caustic network in these images, creating a 
larger difference between the rate of microlensing for sources 
moving with and against the shear. The variation with ori- 
entation of the measurements of effective transverse velocity 
is slightly larger when the model contains a component of 
smooth matter due to additional stretching of the caustic 
network. 

A similar argument applies to the determination of the 
equivalent transverse velocity of a stellar velocity dispersion. 
The rate of microlensing due to stellar proper motions at a 
series of fixed source points is independent of the direction 
in which those points are sampled, regardless of the con- 
tributing fraction of smoothly distributed matter. However 
microlensing that results from a transverse velocity has a 
rate that varies more with direction in images C and D. The 
equivalent transverse velocity is therefore a function of direc- 
tion. The effect is more pronounced in models that include 
a smoothly distributed mass component. 

Thus there are two separate effects. Firstly, if the source 
trajectory is aligned with the image C-D axis the measured 



The calculations in Sec 4.2 and Sec 4.4 assumed that mi 



crolensing results from the combination of stellar proper 
motions and galactic transverse velocity. However the bulge 
dynamics may be dominated by rotational motion result- 
ing in stream motions of the starfield. For an isothermal 
sphere model of the bulge, the rotational velocity is \/2o"* = 
233 fern sec" \ Kundic, Witt & Chang (1993) obtained the 
expression for the caustic velocity Vcaust resulting from a 

stream motion V stream — {Vatream.,^)'- 



^caust 



■ 7 - 



\Vstr 



(10) 



If the rotational motion is assumed to be circular then the 
stream motions are perpendicular to the image-galactic cen- 
tre axis, and therefore parallel to the shear vector at each 
image. We choose the stream motions to be in the xi direc- 
tion making the shear positive for all 4 images. An equivalent 
transverse velocity was defined for the stream motion and 
determined in analogy to that for the stellar proper motions. 
For a circular stream motion of v stream = 233 fcm see" ^ 
equivalent transverse velocities {Vequir{STREAM)) were 
calculated for the two transverse directions discussed pre- 
viously. These values are presented in Tab. ^ When com- 
bined with the equivalent velocities for an isotropic velocity 
dispersion, the values describe the contribution to the mi- 
crolensing rate of the extreme cases for stellar motions in 
the galactic bulge. 

The lower panel in Tab. ^ shows results for models 
with stream motions. Values are shown for mjom (99%), 
miom{95%), miom{90%) and the mode for all models con- 
sidered, mjo^ (99%) is > O.OOSMq in all models. The equiv- 
alent transverse velocities presented in Tab. ^ show that an 
orbital stream motion produces a much lower microlensing 
rate than a transverse velocity of equal magnitude (indepen- 
dent of its direction). The reason for this behaviour is ap- 
parent from Eqn. |l^; a positive shear reduces the caustic ve- 
locity resulting from stream motions. In addition, the caus- 
tic velocity and hence microlensing rate are also decreased 
when there is a component of smooth matter. The lower 
microlensing rate translates to an equivalent transverse ve- 
locity that is smaller than the one obtained from a stellar 
velocity dispersion. This in turn means a lower estimate of 
miom{P)- In the case discussed, the mass limit obtained is 
only ~ 0.2 times that where the proper motions are assumed 
to be isotropic. 

The last row of Tab. ^ summarises the range of values 
for the mean microlens mass obtained by this study. These 
compare favourably with the results of related studies, none 
of which find that Jupiter mass compact objects are the 
dominant component of mass in the bulges and or halos of 
the spiral galaxies studied. 
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5 CONCLUSION 

Calculations of the transverse velocity from published mon- 
itoring data of Q2237-I-0305, as well as complimentary char- 
acterisations of caustic clustering from previous studies show 
that the caustic networks produced by populations of mi- 
crolenses have a structure independent of the mass spec- 
trum, but a scale length proportional to the square root 
of the mean microlens mass. Previous determinations of 
average mass from microlensing data have therefore been 
proportional to the square of the unknown transverse ve- 
locity. However a minimum rate of microlensing is pro- 
duced by the proper motions of microlenses in the bulge 
of the lensing galaxy. We have used limits on the effec- 
tive transverse velocity obtained from published monitoring 
data for Q2237-I-0305 in combination with calculations of 
the microlensing effect of isotropic stellar proper motions 
to calculate the probability for the mean microlens mass. 
We find that the most likely value for the mean microlens 
mass is O.O4-O.45M0. The lower limits on the mean mi- 
crolens mass responsible for the observed microlensing are 
(m) > 0.02 - O.24M0 (99% level). A similar calculation as- 
suming an orbital stream motion of stars rather than a veloc- 
ity dispersion produces a most likely value of O.Ol-O.O5Af0 
and a lower limit of (m) > 0.005 - O.OTIMq (99% level). 
Through joint consideration of the probability for mean mi- 
crolens mass and transverse velocity under the assumption 
of an isotropic velocity dispersion we obtain a mass between 
O.OIMq and I.OMq (95%), and a galactic transverse velocity 
less than 300 fcmsec'^ (95%). 

Our result that very low mass objects do not 
comprise a significant mass fraction of microlenses in 
Q2237+0305 supports the results of Schmidt & Wambs- 
ganss (1998) who rule out microlenses in Q0957-I-561 having 



(m) < 0.01 



- — ^) . The most likely mean mi- 

a sec / 

crolens mass of ~ O.2M0 supports both the conclusions 
of Lewis & Irwin (1996) who found that microlenses in 

Q2237+0305 have O.IMq ^ (m) ( gpo ^ IOMq, 

and of the MACHO collaboration (Alcock et al. 1997a,b) 
who find from microlensing in the Milky Way that (m) ^ 
O.O5M0 in the halo and O.IMq <; (m) <, I.QMq towards the 
bulge. 

While our limits depend on the correctness of the pub- 
lished macrolensing parameters, as well as on the fraction 
of smoothly distributed matter assumed, they do not de- 
pend on other unknowns present in the problem including 
the source size, source intensity profile, effective transverse 
velocity and direction of the source trajectory. We have in- 
vestigated models having smooth matter fractions of and 
0.5, and find that smaller masses can explain the observed 
microlensing rate if a non-zero smooth matter fraction is 
assumed. Our results depend on the Bayesian prior chosen 
for galactic transverse velocity. However we have chosen pri- 
ors to bracket physically reasonable choices and find that 
the contribution of the prior to the systematic uncertainty 
is smaller than that of the assumption for smooth matter 
fraction. 

Our results unambiguously rule out a significant con- 
tribution of Jupiter mass compact objects to the mass dis- 
tribution of the bulge whether the microlenses move pre- 
dominantly in orbital or random motions. However if the 
microlens proper motions are distributed according to an 



isotropic velocity dispersion then the lower limit obtained is 
of the same order as the hydrogen burning limit. 
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